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Page  1,  Abstract,  line  2,  insert  (see  Box  & Jenkins,  1976)  after  "well  known. 
Page  2,  line  4,  E(z  ) should  be  E[z  ]. 

X f u X / u 

Page  2,  Eq.  (1.2),  brackets  are  missing  in  Eta  ^ 

X f u X • n f 

Page  3.  The  following  should  appear  as  a caption  below  the  figure; 


= y y y ii* 

,x_,t  ^ ^ . ,^,’^n,,n_,) 


. a . . , ■'■a  .. 

, ,..„,k  x,+n,  ,x_+n- ,t-k  x,  ,x_,t 

n =-l  nj^=-2  k>=l  12  1122  12 


spatial  order  four  in  direction  x^^,  order  two  in 
"Whether  or..."  should  begin  a new  paragraph. 

Page  3,  line  2 from  bottom  of  page  "at"  should  read  "as." 

n -n^^ 

Page  4,  line  1,  a comma  should  precede  F and  a comma  should  follow  B 

X X- 

m m 1 

Page  4,  line  14,  should  read  X S^, 

i=-m 

2 2 

Page  9,  insert  following  Eq.  (2.1.9)  "Further  since  ^^lO^a^ '^'^z' 

then  we  find  terms  of  and  only. 


-Pll-Pll 


yi-4(p2,+p2,) 


01  ^11 


(2.1.10) 


The  advantage  of  (2.1.9)  over  (2.1.10)  is  that  6^^  and  are  unambiguously 


determined  by  p^^^,  p^^^,  and  p^^^. 


Also  change  original  (2.1.10)  to  (2.1.11). 


Page  10,  following  (2.2.3)  insert  "Note  P_j^j^~Pj^]^'  because  of  the  sym- 
metry of  the  autocorrelation  function  in  m dimensions." 

Page  11,  line  15,  add  after  "implies"  "|6^|<1,  again  in  agreement. 
Finally,  with  obtain  the  MA(2;0,0)  MA(2)  process  of  Box 

and  Jenkins,  and  (2.2.5)  implies..." 

Page  11,  line  5 from  bottom,  add  "In  (2.2.3)  Pj^2.**^-ll  symmetry 

of  the  autocorrelation  function  in  m dimensions." 
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"Note 

^llo"^l-10 

Page 
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line  12,  after  "zero" 
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"Note 

' 101 '^-101'  ^011  ^0-11' 

^110=^1-10- " 

Page  20,  line  13,  following  "autocorrelations"  add  "Note  that 


'^lor'-ioi'  'oirVii'  '^iio  Vio'  ’^iir''-iii'  ^i-ir’-i-ir 
Page  21,  line  17,  should  read  "eralize  to  m dimensions  under  proper 
restrictions. " 
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Taneja,  V.,  & Aroian,  L.A.  (1977).  Time  series  in  m dimensions, 
autoregressive  models. 
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ABSTRACT 

Stochastic  models  for  discrete  time  series  in  the  time  domain 
are  well  known  but  such  models  lack  consideration  of  spatial  de- 
IJendency.  We  expand  on  their  work  by  constructing  spatially  depen- 
dent moving  average  models.  Definitions  of  order,  stationarity , 
invertibi lity , autocorrelation  function,  and  spectrum  are  made  as 
natural  extensions  of  those  in  zero  dimensions  and  are  implemented 
in  the  one  and  two-space  dimensional  models. 


1.  INTRODUCTION 

We  describe  a general  linear  stochastic  model  which  supposes 

a time  series  to  be  generated  by  a linear  aggregation  of  random 

shocks  at  various  temporal  and  spatial  locations.  Letting  x = WOH  for 

(x,,x_,...,x  ),  an  m-dimensional  vector,  the  general  Moving 
12  m 


Average  (MA)  model  of  m-dimensional  time  series  is  defined  by 


Wjijia  Stcflwi 
F-.i?  Ssctlon  □ 


JUSil.  I'll'  { 
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^Xjt  ^ ^ j^^j^'^n,k*x+n,t-k  ^ ®x,t' 


where 


n = (n, ,n„,...,n  ),  J denotes  the  repetiti 
12  m _ 


ive  sum  over  each 


00  00 


component  ofn(i.e.,  'l  = ^ •••  I I ^ = z - 

X X #c 

n=-oo  n =-"  n^=-«>  n =-" 

m 2 1 

E(z  ) is  the  deviation  from  the  mean.  The  white  noise  process 

X rt 

a may  be  regarded  as  a series  of  independent  random  shocks  which 

X /t 

drive  the  system,  so  that  their  autocovariance  function  is: 


Y ,=Ea  .a^  .. 

n,k  . x,t  x+n,t-k 


0 , n=0  and  k=0 

a 


0,  otherwise 


and  hence  the  autocorrelation  function  of  white  noise  has  the  par- 


ticularly simple  form 


1,  n=0  and  k=0 


n,k  10,  otherwise 


We  also  assume  that  z is  a weakly  stationary  process,  i.e., 

X ^ t 

Etz^  ] < a.  and  e[z  . z.  . ] =»  a ^P|  11  I . (1.4) 


x,t^  y.t^ 


|x-y| , 


In  this  paper  we  explicitly  focus  our  attention  on  the  special 
case  of  (1.1)  in  which  only  a finite  number  of  the  coefficients  are 


nonzero,  that  is: 


*x,t  j^^j^'^n,k*x+n, t-k  ^ *x,t 


If  in  equation  (1.5)  all  of  the  coefficients,  ip  , are  nonzero, 

n f K 

tlie  process  is  called  moving  average  of  temporal  order  r and  spa- 
tial order  p^  + in  each  space  direction  x^,  1 < j < m.  For 
example,  letting  m = 2,  p = (2,1),  q - (2,1),  and  r arbitrary,  we 
have  the  following  representative  scheme; 
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!.paticjJ  ortifr  four  in  direction  x^,  order  two  in  x^.  Whether  or 
not  all  the  coefficients  are  nonzero,  it  is  easier,  from  the  anal- 
ysis viewpoint,  to  represent  the  process  (1.5)  in  terms  of  shift 
operators.  The  backward  shift  operator  in  time,  B^,  is  defined 
by 

®t^x,t  “ ^x,t-l' 

while  the  backward  and  forward  shift  operators  in  spatial  direc- 
tion X.,  denoted  by  B and  F respectively,  are  defined  by 

^i  i 

B z — z p 
X.  X, t x-6. ,t 
1 1 


(l<i<M) 


F z = z , , 

x^  x,t  x+0^,t 


where  6.  = (6.,,  6. .,•••»  6,  )and  *) . . = 
1 il  x2  im  13 


i=j 

0, 


Powers  of  these  operators  are  defined  in  the  usual  manner,  for 
exeunple, 

2 - - _ 

B z ^ = B (B  z .)  = z _P 
X.  x,t  X.  X.  x,t  x-26.,t 

111  1 


In  addition,  wo  note  that  the  operator  B is  t)ie  inverse  of  F , 
that  is,  B^  = . In  terms  of  these  shift  operators,  the 


process  (1.^)  can^be  reformulated  at 


= (1+  ^ I 'i'  B*')  - 

c,t  n=-p  k=l  ^ ^ 


) a 


(1.6) 
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Defining 


^(B  ,B  I = I * f J 
n=-p  k=l 


the  moving  average  process 


(1.7) 


z'  ^ = '('(B  ,B^)a  ^ 

X , t X t X , t 


(1.8) 


can  Ix!  thought  of  as  the  output's  from  a linear  filter  with 

X / t 

transfer  function  y(B  ,B  ) when  the  input  is  white  noise  a 

X t x,t 

Since  the  expression  for  'l'(B^,B^)  is  finite,  no  restrictions  are 
needed  on  the  parameters  ip  to  ensure  stationarity . The  in- 

n f K 

vcrtibility  condition  for  the  moving  average  process  may  be  ob- 
tained by  writing  (1.8)  as 

a ^ = 4'‘^B  ,B^)z 
x,t  x t x,t 

Extending  the  results  of  Box  and  Jenkins  it  can  be  shown  that,  for 

-1 

invertibility , II(B  ,B  ) =4'  (B  ,B^)  must  converge  on  X S^, 
where 


= {B^:|bJ<1} 


s. 

1 


(l$i<m) 


The  autocovariance  function  of  a MA  process  may  be  obtained 

by  multiplying  through  (1.6)  by  “z  . , , where  ii  = (H  , Z , , i ) 

J.  A in 

and  talcing  expectations.  A more  convenient  way  of  obtaining  the 
autocovariances  is  often  via  the  autocovariance  generating  function 


■ I I 

00  GO 


(1.9) 


in  which  it  is  noted  that  Yqq»  the  variance  of  the  process,  is  the 
coefficient  of  B^B*^  z 1 , while  y„  , is  the  coefficieit  of  both 

X t ?,,k 
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I 
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! 

i 

f 
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1 

1 


I 


i 
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p 

B 13  and  of  B B ; as  in  Box  and  Jenkins,  it  follows  that 

X t X t 

Y(B  ,B  ) = o'  4'(B  ,B  ) I*  (F  ,F  ) . (1.10) 

X t a X t X t 

The  autocorrelation  function  in  m dimensions  is  symmetrical  in 
each  variable. 

-i2TTf 

It  can  be  shown  that  if  we  substitute  B^  = e and 

B = o whore  i = and  l.^j^m,  in  the  autocovariance 

X . 

3 

generating  function  (1.10),  we  obtain  the  power  spec- 
trum. Thus  the  spectrum  of  the  MA  process  is 


p(g,f)=2a 


2 

a 


,,, , -i2TTg  -i2'rrf  i2Trg  i27Tf 
f(e  ’,e  )4'(e  ’,e 


(1.11) 


-i2nf  I 2 

e ) I , 


where  0<fiS,  0Sg^.<S  (l5Cjim)  , and  4'(e  ^^^'^)='('(e  ^^^^1 , 

-i2'nqp  -i2TTgn,  -i2'nf, 

e e , e ) . 


In  sections  two  and  three  we  analyze  one-dimcnsional  and  two- 
dimensional  models,  respectively.  In  both  sections  we  have  fo- 
cused on  models  whose  spatial  and  temporal  orders  do  not  exceed 
two.  Section  four  contains  some  general  remar)cs  and  indicates 
areas  where  more  extensive  research  is  needed.  The  following  the- 
orem generalizes  a corresponding  result  for  MA  in  zero  dimensions. 

Theorem:  If  the  conditions  for  invertibility  are  satisfied,  then 

every  finite  MA  process  in  m dimensions  may  be  ex- 
pressed as  an  infinite  AR  model  in  m dimensions. 

2.  ONE-DIMENSIONAL  TIME  SERIES 

In  this  section  we  consider  moving  average  models  of  one- 
dimensional time  seri« '■  with  temporal  and  spatial  orders  not 
exceeding  two.  foj  .'nvenience  and  easy  reference  we  let  x=Xj^, 
and  denote  by  MA  (r;Pj^,qj^)  the  moving  average  model  for  one- 
dimensional  time  series  of  temporal  order  r and  spatial  order 
Pj^+q^;  recall  that  p^  and  q^^  denote  the  maximum  powers  of  the 


6 


operators  and  occurring  in  (1.6),  respectively.  This  nota- 
tion is  an  attempt  to  be  consistent  with  that  in  Box  and  Jen)cins, 
who  analyzed  zero-dimensional  time  series  with  moving  average 
models  of  various  temporal  orders,  that  is,  MA (r )~MA (r ; 0, 0) ; this 
equivalence  means  that  since  the  maximum  powers  of  the  spatial 
operators  are  both  zero,  and  B^  - = 1 (identity  operator),  we 

are  virtually  loo)cing  at  the  same  point  at  different  times. 
Formally,  from  (1.6)  we  find  that 


"x,t 


r 

I ^ 

n^=-p^  k=l 


(1+  I 


*^1  k 

n^,)<^x  ^t^  ^x,t‘ 


(2.1) 


Thus  with  Pj^=qj^=0,  we  obtain 


^x,t  = ‘*x,t^^l  ^,t-2^-"^V  ‘^x,t-r- 


Deleting  the  first  subscript  in  each  case,  as  z depends  only  on 

X # c 

t,  results  in  the  zero-dimensional  moving  average  model  in  Box  an<i 
JenJtins,  namely; 


z 

t 


a +il;  a + a +. 
t ^1  t-1  2 t-2 


.-t-ip  a 
r t-r 


(2.3) 


There  are  five  general  "spatial"  models  described  by  (2.1) 
where  lSPj^+qj^<2;  we  represent  these  diagramatically  below,  where 
order  refers  to  the  spatial  order: 

x-1  X 


MA(r;l,0) : 
MA  (r ; 0, 1)  : 
MA(r;l,l)  : 
MA(r;2,0) : 
MA(r;0,2)  : 


First  order  backward 
First  order  forward 
Second  order  forward-backward 
Second  order  backward 
Second  order  forward 


X 

x+1 

x-l 

X 

X4  1 

x-2 

x-l 

X 

X 

xn 

x+;‘ 

As  indicated  previously,  we  will  restrict  our  attention  to  the 
cases  r=l,2;  this  results  in  a set  of  ten  models  of  which  we  will 
analyze  a subset. 
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2.1  The  model  MA(1;1,0).  From  equati'  ■ (2.1)  we  obtain 

0 1 n 

B ) a 

,t 


* * * 1 

s *.=  (1+  I I 'I'  uF  B^)  a 

x,t  ^ n^,k  X t X, 


(2.1.1) 


= a +ili  a .+ii)  a 

x,t  ^01  x,t-l  ^-11  x-l,t-l 

For  convenience,  and  to  indicate  that  we  are  using  a finite  set 
of  weight  parcuneters,  we  change  symbols,  letting  and 

~”®2‘  (2.1.1)  becomes 

z ^=(1-O,B^-0^B  B^)a  ^=G(B  ,B^)a  (2.1.2) 

x,t  It  2xt  x,t  X t x,t 

Multiplying  through  (2.1.2)  by  z o ^ 

X”  X/  f K 

^x-£.,  t-)c^x,  t ^x-)!,,t-)<^x,  t ^l^x-?,,t-)c^x,t-l  ^2^x- £,  t-k^x-1 , t-1 

(2.1.3) 

On  taking  expected  values  in  (2.1.3)  we  find  the  variance  of 
the  process  is 

wv 

and 

^01=-®l”f'  ^10=°l®2^^  ^ll=-'^2‘'f' 

while  all  other  autocovariancos  are  zero.  Thus  the  autocorrelation 
function  is 

-01  0^02  -02 

^Ol"  '2  2"*  Pio"  ~ 2 2 ' Pll"  r~2'  <2.1.5) 

1+01+02  ^^^^^2 

all  other  being  zero. 

To  illustrate  the  technique  of  obtaining  the  autocovariances 
using  the  autocovariance  generating  function  (1.9), 


y(B  B^)=0(B  ,B^)=l-0  B^-0.B  B^ 
xt  xt  lt2xt 


(2.1.6) 


I 
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is  suhsl  i t ul  i'(i  into  (I.IO)  yicldinq 

Y(B  ,B^)=O^(1-0,B^-O  B B^)  (1-0  (2.1.7) 

xt  a lt2xt  It  2xt 

= 0^  [-0  b'^b'^+0,  0,b"^-O  b"^+(1  + 0J+0^) 

a 2xt  12x  It  12 

-0  B +0  0 B -0  B B ] . 

It  12x  2xt 

Comparing  (2.1.7)  witli  (1.9)  and  noting  that  Y_i-i  “Yu, 

Y ~^10'  ^0-1  ~''*"oi'  obtain  the  same  results  as  listed  in 

(2.1.4) . 

Substituting  i2Trg  B^=e  into  the  autocovariance  gen- 

erating function  (1.9),  with  ^(B  ,B  ) as  in  (2.1.6),  we  obtain  one 

X t 

half  of  the  power  spectrum.  Thus  the  spectrum  is 


p(g,f)=2a^|4'(e’^^’"^,e‘'^^’'^)  (2.1.8) 


_ 2|,  a “l2Trf  a -l2-ng  -l27Tf  i2 
=20  1-0, e -e_e 

a 1 2 ' 

= 2o^  ll+0^+6a+2O  0 cos27ig 
d X ^ X ^ 


-2  {0^cos2Trf+02COs27r(f+g)  }] 
0<fS*3;  0igl*5 


For  invertibility , the  generating  function  n(B^,B^)=0  (B^,B^) 

must  converge  for  |b^|<1  and  |B^|i^l.  Since 

OO 

3=0 

wo  see  that  for  invertibility, 


|01+O2Bxhl  and  |b^|.<1. 


Consequently,  the  parameters  of  the  MA(1;1,0)  process  must  satisfy 
|0ll+|02l<l  to  ensure  invertibility.  If  the  autocorrelations  are 
)«nown  we  can  solve  for  the  parameters  0^^  and  6^  from  equation 
(2.1,5)  even  though  they  are  nonlinear;  for  other  models  we  may 
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need  to  use  an  iterative  process.  In  the  present  case. 


0 


1 


(2.1.9) 


If  any  of  the  autocorrelations  are  zero,  then  from  (2.1.5),  0^ 
or  0^  is  zero,  the  model  reduces  to  the  zero  dimensional  case 
and  we  need  to  solve  a quadratic  equation  for  the  nonzero 
parameters. 

Finally,  we  wish  to  point  out  that  the  analysis  of  the  forward 
model  MA(1;0,1)  given  by 


x,t 


(1+ 


1 1 

I 

=0  )c=l 


"l  k 

,F  ^B^)a  . 

nj^,)c  X t x,t 


(2.1.10) 


=a  +11^  a +i0  a 

x,t  ^01  x,t-l  ^11  x+l,t-l 


is  very  similar  to  that  of  the  MA(1;1,0)  model.  The  results  are 
analogous  with  the  par£unter  0^  replaced  by  0^  and  the  operator 
replaced  by  F^. 


2.2  The  model  MA(2;1,0).  From  (1.6)  we  obtain 


1 V 

? =(1+  y l B^)a  ^ 

x,t  n,  ,)c  X t x,t 

nj^=-l  k=l  1 


(2.2.1) 


=\,t^'<'01  ’x,t-l*'<'-llVl,t-l^''^-12Vl,t-2^'<'02%,t-2 

Where  — 0^,  =-02-  =-6^,  and  =-d^. 
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Substituting  into  the  autocovariance  generating  function  results 


Multiplying  out  this  expression  and  comparing  the  result  with  (1.9) 
we  find  the  variance  of  the  process  is 


Y =O^=(l+0^+9^+0^+0^)O^ 

00  z 1 2 3 4 a 


(2.2.2) 


"01=<-®l^°2V°l'4>"a''  \0=^°lV°3°4’'^f' 
Yii=(-e2te^03)a2,  Y.,,=02'4"^ 

while  all  other  autocovariances  are  zero.  Thus  the  autocorrelation 
function  is 


1.02.02^02.02  1,02^0^.02.02 


(2.2.3) 


-0  .0  0 

-^-4  3 — 

1.0%o2.02+02 

12  3 4 


1.02.O2.O2.O2 


_ -»4  , -O3 

1+0^  i+oJ+gJ+gJ+o^ 

i2j4  1234 


From  (1.11),  the  power  spectrum  for  this  model  is  given  by 


/ rv  - -i2Tlf  „ -i2TTg  -i2TIf  ^ -i2Trg  -i4liC  „ -i4nf|2 

p (g, f )=20^1 l-03e  -d^e  ’e  -G^e  ^e  -0^e  | 


0<f<S;  0<g<S 


(2.2.4) 


il 


For  invertibility,  11(0^ , B^)  =0  ^ (B^ ,B^)  must  converge  for 
|b^  I <1  and  |b^  1^1.  Equivalently,  the  roots  of 

must  be  outside  the  region  |B^|il  and  |b^|^1.  Hence  the  param- 
eters must  satisfy: 


e.-e -0 +e.si,  -e, +9,-0 +0.^1, 

1234  1234 


(2.2.5) 


103  1+104  1^1- 


Note  that  as  a consequence  of  these  inequalities  it  follows  that 

|0J+|02  1^2. 

Setting  02“04~O»  obtain  the  model  MA(1;1,0)  discussed  pre- 
viously. From  (2.2.5)  with  03“04“0  obtain  1 0 1 * I 0^  I 1 , wliii.h 
agrees  with  the  invertibility  condition  for  the  MA(1;1,0)  model. 

In  addition,  note  that  with  02“03”04“0  obtain  the  MA(1;0,()) 
~MA(1)  process  of  Box  and  Jenlcins,  and  (2.2.5)  implies 

01+ 04-01^1'  |04l^l' 

which  is  also  in  agreement.  If  in  (2.2.3)  the  autocorrelations  are 
Itnown,  we  get  a nonlinear  system  of  six  equations  to  solve  for  the 
parameter  0^,  l^ii4.  Certain  dependencies  exist  among  the  auto- 
correlations as  in  the  previous  model,  and  the  resultant  system 
may  be  solved  iteratively. 


2.3  The  model  MA(1;2,0).  From  (1.6)  we  obtain 

0 1 n . 

z ^=(1+  I J ^BJa  ^ 

x,t  , ,n,,)<x  t x,t 

n,=-2  )i=l  1 


(2.3.1) 


'’^x,t'^'^01  ®x,t-l‘*’'^-ll  ®x-l,t-l'*'^’-;M  ‘‘x-2,t-J 


^x,t  ^l^x,t-l  ^2^x-l,t-l  ^S‘’x-2,t-l 
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-0,B  B -0-B  B^)a  ^ 
It  2xt  3xt  x,t 


where  ~”®l'  *^-11  “”®2'  *^-21  ~~®3'  Substituting  into  the 

autocovariance  generating  function  we  find 


Y(B  ,B  )=a^(l-0  B.-0.B  B^-0-B^B^) 
xt  a It2xt3xt 


Multiplying  out  this  expression  and  comparing  the  result  with  (1.9), 
we  find  the  variance  of  the  process  is 


and 


2 2 2 2 2 


(2. 3.2) 


^20'°1^3^a' 


while  all  other  autocovariances  are  zero.  Thus  the  autocorrela- 
tions are  given  by 


-0. 


^01° 


“ I+0J+O2+03 


(2.3.3) 


«lV»2«3 


-0, 


2 ' ^21^ 


2 2 2 
l+0^+0%03 


P 


20 


^«3 


l+0ft0^+0 

1 2 


2 

3 
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p(g 


From  (1.11),  the  power  spectrum  for  this  model  is  given  by 

(2.3.4) 


O^fS*!;  OigcS- 


For  invertibility,  FI  (B^,B^) =0”^ (B^, B^)  must  converge  for 
|b^|s1  and  |b^|^1.  Since  Q(B^,B^)=1- (0^+O2B^+03B^)B^,  inverti- 
bility is  ensured  if  jO  +0,B  +0  B^|<1  and  |b  |<1.  Thus  the  param- 

X ^ X X X 

eters  must  satisfy 

|0i+O2+Q3l<i4  |o^-02+e3l<i. 

0 02 
and  if  — < 1,  then  |0^-  -^1  <1. 

2.4  The  model  MA(1;1,1).  From  1.6  we  obtain 


n,“-l  K=i  i 


(2.4.1) 


"^x.t^’^Ol  *x-l,t-l^'^ll  ^x+l,t-l 

=“x.t-®i%,t-r^Vi,t-r®3Vi,t-i 

■'i-Vt-W4-“3^‘"t’s,t 


where  '*'-11  ^~^2’  *''**  **'ll  “~^3‘  Substituting  into  the 

autocovariance  generating  function 

,(B,,B^).C.^l-eiB^-0jBxB,-e3B;‘B^) , 

and  comparing  with  (1.9),  we  find  the  variance  of  the  process  is 


and 
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>01=- VI' 


\i=-®2^;' 


\o=^®iV®iV°f'  Y-ii=-Vf' 


Y =9  0,0  , 
20  23a' 


(2.4.2) 


while  all  other  autocovariances  are  zero.  Thus  the  autocorrela- 
tion function  is 


-9. 


-0- 


^01 


P,^= 


On- 


'■-n- 


-0, 


i+e^+e^+o^  ' 


®2S 


20  i.o2.e2.02 


From  (1.11),  the  power  spectrum  is  given  by 


(2.4.3) 


, ,_2  I,  n -i2Hf 

p(g,f)=20^  |l-0^e 


-i27Tq  -i2Tlf  f,  i2i\q  -i27rf  |2 
■O^e  "'e  -O^c  e | , 


(2.4.4) 


0<:f<‘j;  OSgSb. 

.1  1 

For  invertibility , n(B^,B^)=0  (B^,B^)  must  converge  on  X S^; 

since 

0(B  ,B^)=1- (0  +0_B  +0,F  )B^,  we  obtain 
X t 1 2 X 3 X t 

101+02+03  I<1»  |0i+O2-®3l^^' 

101-02+03 I<1.  101-02-03  I<1. 

or,  equivalently,  |0j  |+ [02  1+ I03  Hi  conditions  on  the  parampt<;rs. 

3.  Two-Dimensional  Time  Series.  In  this  5:oction  we  consider 
MA  models  of  two-dimensional  time  scries  with  tt*m|>oral  and  sp.itial 
orders  not  exceeding  two.  For  convenience,  we  let  x«x^,  y“yj  and 


I'j 


denotu  by  MA(r;p,q),  where  P=(p^/P2)»  ^ process 

of  temporal  order  r and  spatial  order  x-direction 

and  V>2*^2  y-direction. 

From  (1.1)  the  general  model  has  the  form 

OO  00  00 

~z  .=  y y y ij;  i,® . .4.  o.d 

x,y,t  ^ ,^,^n  ,n_,k  x+n,  ,y+n  ,t~k  x,y,t 

n_=-“’  n,=-‘”  k=l  12  12 

2 1 

The  special  case  of  (3.1)  in  which  only  a finite  number  of  the 
coefficients  are  nonzero  results  in 

q,  qi 

^'x,y,t  ^ , ^.'^n,  ,n_,k^x+n  ,y+n  ,t-k^^x,y,t' 

n^=-P2  12  12 

or  in  terms  of  shift  operators, 


^*2  ‘^1 


"l„"2„k. 


= (1+  y y J rli  F F B )a  . 

,y,t  '■  ^ , n,  ,n_,k  x y t x,y,t 

n„=-p^  n =-p,  k=l  1 2 


n2=-P2  n^=-p^ 

= '1'(B  ,B  ,B  )a  , . . 
X y t x,y,t 


In  addition  to  the  zero-dimensional  model  corresponding  to 
p=q=0,  and  the  ten  one-dimensional  models  in  each  direction  re- 
sulting by  letting  Pj^=qj^=0  or  there  are  fifty  general 

"spatial"  models  described  by  (3.3)  where  l<p^+q^';2,  i=l,2,  and 
l$r-^2.  In  what  follows,  we  examine  three  such  models. 


3.1  The  model  MA(1;1,1,0,0) . From  (3.3)  we  obtain 


0 0 1 


"l  ”2„k. 


-.y.t'"*  ^ ^ V\,y,t  (3.1 

n^='”l  K=1  1 2 

"®x,y,t^'^  001  *x,y,t-l''^'^-101  ®x-l,y,t-l 

■^'^0-11  ®x,y-l,t'l'^’*'-l-ll  ®x-l,y-l,t-l 

“*x,y,t  *^l®x,y,t-l  ®2®x-l,y,t-l  ^3®x,y-l,t-l 

“S®x-l,y-l,t-l 


= (1-0,B^-0.,B  B^-O.B  B^-9^B  B B ) a 

It  2xt  3yt  4xyt  x,y,t 


i 
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where  tooj  -6^,  -Sy  Vll  »'''“  ♦-1-11 

Substituting  into  the  autocovariance  generating  function 

Y(B^ . By , Bt ) . l-  e^B^-  0^B__B^-  638^6^-  r 

(1- 9, b"^- 0^b"^- 0,b”^b*^- 0 b"^b"^b"^) 

It  2t  3yt  4xyt 

and  comparing  with  (1.9),  we  find  the  variance  of  the  process  is 


and 


''ou■-V^  ''no*  '’iVe 


(3.1.2) 


while  all  other  autocovariances  are  zero.  Thus  the  autocorrf 
lation  function  is  given  by 


-0. 


«1«2  ■^®3«4 


001 


101 


oil 


-0- 


1 + 0S0S0S0^ 

12  3 4 


-0, 


-0. 


l+0j+o5+0^+0f 

12  3 4 


' ^100' 

2222' 

1^0^02^03^+0,' 

^V®2®4 

' ^010° 

1+0^+0^02^02  ' 

®1®4 

' ^110 

«2«3 

' ^1-10' 

i+6i+o^G^o; 

(3.1.3) 


-i2TTgi  -121192  -i2ltf 

On  substituting  B^=e  , By=e  , and  B^=e  into 


(1.11),  we  find  tlie  power  spectrum  of  the  process  is  given  by 

(3.1.4) 


2,  -i2iTgi  -i2iTg 

p(g,f)=2cJ  l-(0  -0,e  ^-0  e .. 

a 1 2 3 4 


-i2iigi  -i2lig2  -i2iif|2 
-0  e e )e 
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For  invertibility , 


ll(B  ,B  ,B^)=0  (B  ,B  ,B^) 
X y t X y t 


= [l-(9  +0^B  +0.B  +0.B  B )B  ] 
1 2x  3y  4xy  t 


-1 


miiRt  c.'onvoryc  on  X S..  Hence,  the  parameters  must  satisfy 
i=-2  ^ 

>VW®4l'^'  lVW®4l^^' 

to  ensure  invertibility. 

3.2  The  model  MA( 1 ; 1, 1, 1 , 1) . From  (3.3)  we  find 


111 
.=  (1+  I I I 


"l  "2  )c. 

F F B ) cl 

'x,y,t  '■  , v-1 X y t x,y,t 

n2==-l  n^=*i  K-i  i 2. 


■^'^-ll^x+l,y-l,t-l‘^'^-101®x-l,y,t-l  (3.2.1) 

‘^'*001^x,y,t'*’'^101®x+l,y,t-l'^'^-lll^x-l,y+l,t-l 

^'*'oil®x,y+l,t-l'^'^lll®x+l,y+l,t-l’ 

Consider  the  special  case  where 

Letting  ^01=-®3'  Vll=-®4  ^ll=-®5' 

we  obtain 

^x,y,t  ^x,y,t  *^l*x,y,t-l  2^x-l,y,t-l  3 x+l,y,t-l 

^^4*x,y-l  ,t-l  *^5^x,y+l  ,t-l 

= [l-(0  >0.,B  +0,b"^+0.B  +0.B"^)B^la 

12x3x  4y5y  t x^y^t 

=0(U  ,B  ,B^)a 

X y t x,y,t 
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Substituting  into  the  autocovariance  function 

Y(B  ,B  ,B  )=a^  (1-0,B^-0^B  B^-0,B  "^B^-0  B B^-0_b"^B  B)x 
* x^y^t  d It  2xt  -ix  t 4yt  ^y  yt 

(1-0  b“^-0^b"^b"^-0,b  b'^-0^b"^b"^-0^b  B-1, 
It  2xt  3xt  4yt  5yt) 

and  comparing  with  (1.9),  we  find  the  variance  of  the  process  is 


L=-®i'’2 

W>l<V«3’'’f 

L=-®2"a 

=-Va 

^200-'’2’‘A 

II 

1 

‘^020“’4’5’a 

=-0,0^ 

5 a 

<U0*'*2V“3'’4»’; 

II 

c 

1 

(6204*0395)0^, 

while  all  other  autocovariances  are  zero.  From  these,  we  can  ob- 
tain the  autocorrelations. 

From  (1.11),  the  power  spectrum  is  given  by 

, -i2Trgi  ^ i2'ng,  -i2Tkj-,  i2-ng2. 


p(g,f)=  ll-O^  + e^e  +03e 


(3.2.3) 


For  invertibility , 

ri(B  ,B  ,B^)=Q“^  (B  ,B  ,B^) 
X y t X y t 


Il-(0  +0_B  +0  F +0.B  +0  F )B^] 
1 2x  3x  4y  5y  t 


must  converge  on  X S . ; hence  we  need 

i-2  1 

l'’l*«2V<'3'’x*^’y*Vyl'=‘ 


|B^|<l3  |F^|<1»  iBylH,  |Fy|<l  . 

Substituting  in  the  values  il  for  the  (dummy)  variables  results 
in  a system  of  inequalities  that  must  be  satisfied  by  the 
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parameters  to  ensure  invertibility . Aroian  finds  the  values  of 
0.  in  terms  of  the  essential  correlations  » P,r,  » ’ 

P„^  , and  o . 

02.  a 

}.3  The  model  MA (2 ; 1, 1, 1, 1) . We  consider  only  a special 
case  of  thi.s  model  consisting  of  the  model  in  section  3.2  except 
with  a lag  of  two  in  time.  We  obtain 
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Y = (-0  +0  0 +0  0 +0  0 +0  0 +0  0,^)0^, 

001  116  27  38  49  5 lo'  a' 

Y =(0  0 +0  0 +0  6 +0  0 )0^,  Y ,^=(0,0  +0,9.  + 0,6„+0„9,o)Cf^ 

100  '121367  68a  1-10  '24  35  79  8 10  a 

Y =(0  6 +0  6 +0  0 +0  0 )a^,  Y, , =(0,0„+0^0_)a^, 

010  ' 1 4 1 5 6 9 6 10^  a'  111  ' 3 9 5 7'  a' 


"-101=<-S^®2'6^®l®8>"f' 

^11=  ‘-^^®l®9^^®6>°a' 

''^002"~®6‘^a' 

^200-‘®2®3*«7«8'‘’f’ 


Y =(0  0 +0  0 )a  , 
-111  '29  5 8a' 


^l-ll"^®3®10'*'®4®7’'’a 
Y_i_ii=  ^02®1O'^®4®8^°< 
^012^”  ^®9'^®10^'^a 
''^02l''^®4®10‘^®5®9^°a 
^02■-'«7*«8’“^ 
■'20l"®2®e*<’3®7’“f' 


(3.3.3) 


all  other  autocovariances  being  zero;  from  these  we  can  determine 
the  autocorrelations. 

From  (1.11)  the  power  spectrum  is  given  by 


p(g,f)=2a2|l-(0j^+02e 


- (0^+0^e 


-i2^gi  i27Tgi  -i2TTg2  i2TTg2  _2Tjf 
+636  +04®  +656  )e 

-i2TTg^  i2TTgi  -i27rg2  i2iTg2 


+0106 


4Ff , 2 


For  invertibility,  II(B  ,B  ,B  )=0''‘'(B  ,B  ,B  ) must  converge  on 
2 X y u X y t 

X S..  Alternatively,  the  roots  of  0(B  ,B  ,B  ) must  be  outside 
i=-2  ^ X y t 


X S.,  and  we  find 
i.-2  ‘ 

V®7=x*Vx*®8“,*®10V®l*®2'x*®3’'x*®4“y*“5^‘^’ 

'>6»«7“x*»8Wy*‘'l0V®r*2‘x-®3V^V®5'y‘'' 


0,+e-B  .6„F  .0^F  +0,„F  I<1. 
6 7x  8x  9y  10  y' 
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r:cl  I iii<|  I h<.*  Viir i otjual  to  *1  Ip.ifl;;  to  a of 

iiK  <4ual  it  ios  in  tlic  ton  p^r^eters . 

4.  Conclusions.  We  have  defined  the  properties  of  a gen- 
eral moving  average  mfidei  of  m-dimensional  time  series.  Specific 
one-dimensional  and  two-dimensional  models  have  been  investigated. 
Further  work  remains  to  be  done  on  moving  average  models  of  n time 
series  in  m-dimensions. 
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